Physical interpretation of the quantum two-stream instability 
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Abstract 

The unexpected features of the two-stream instability in electrostatic quantum plasmas are inter- 
preted in terms of the coupling of approximate fast and slow waves. This is accomplished through 
the factorization of the dispersion relation into different sectors having positive or negative energy. 
Therefore, the concept of negative and positive energy waves is useful not only for classical, but for 
quantum plasmas as well. We discuss the limitations of the quantum two-stream model in view of 
the weak coupling assumption. 



PACS numbers: 52.35.Qz, 41.75.-i, 03.65.-w 



I. INTRODUCTION 



Since it has been discussed for the first time in the framework of the quantum hydrody- 
namical model l|, the quantum two-stream instability has attracted considerable attention 
in the literature. The reason for this is that it is a benchmark displaying many of the parti- 
cularities of quantum plasmas, including a new unstable branch of the dispersion relation for 
large wave-numbers and almost stationary, quasineutral, nonlinear oscillations 1| without 
analog in classical plasmas. Later on, a kinetic (Wigner-Poisson) treatment of the problem 
[2I showed that temperature effects can suppress the quantum instability. The quantum fluid 
equations have been further applied to several quantum streaming instability problems, like 
the three-stream quantum plasmas 3|, the quantum dusty plasmas 4j, the electron-positron- 
ion quantum plasmas or the magnetized multi-stream instability [6|. The hydrodynamic 
brmalism has also been applied to the quantum filamentation instability, with 7| or without 
8| magnetization. 

The purpose of the present contribution is to provide an intuitive explanation of the 
quantum two-stream instability, in terms of the coupling of electrostatic modes with distinct 
energy contents. Indeed, literature still lacks a physical ("with the hands") understanding 
of the quantum streaming instabilities. In addition to analytic and numeric approaches, the 
physical interpretation of quantum plasma effects is a most welcome extension, helping to 
describe the unexpected features of these problems. As will be shown in the sequel, it is 
possible to identify approximate positive and negative energy modes in electrostatic two- 
stream quantum plasmas. Moreover, the energy exchange between these waves provide the 
clue for the stability analysis in such systems. 

In addition, we also examine the range of validity of the quantum two-stream model. As 
shown in Section IV, the quantum parameter H, which is the ratio between the plasmon 
energy and the kinetic energy of each beam, should be small in order to fulfill the collisionless 
assumption. As far as we know, a clear statement about the smallness of H is missing. 
Nevertheless, the more intriguing aspects of the quantum two-stream instability appear for 
small H, where the so-called quantum instability branch for large wave-numbers appear 
The present contribution gives a physical justification for this exotic unstable quantum 
branch in terms of negative energy waves. Finally, the fact that a large quantum parameter 
tends to enhance collisions is a reason more to develop a (still lacking) efficient collisional 
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quantum plasma theory. 

Negative energy modes are a well-known tool for the analysis of streaming instabilities in 
classical systems p.[lol]. In this work we show that the same ideas are also useful to quantum 
plasmas. The basic heuristic concept of negative energy wave is as follows. Consider, for 
definiteness, coherent wave motion taking place in some dispersive medium crossed by a 
number of streaming particles. When the absolute value of the wave's phase velocity is 
slightly smaller than the speed of some of the other streams, the oscillation mode is referred 
to as a negative energy wave. The mechanism of the classical two-stream instability can be 
explained in terms of energy transfer from a positive energy Langmuir wave taking place in 
one beam, to a negative energy Langmuir wave taking place in the second beam. Naturally 
in quantum plasma physics such a classical picture of resonant wave-particle interaction 
should be taken with caution. However, we will show that negative energy modes can still 
be identified in terms of macroscopic (fluid) quantum plasma equations for the two-stream 
instability. 

This work is organized as follow. In Section II, we review the basic features of the 
quantum two-stream instability, and of the electric field energy density in a weakly dissipative 
dielectric medium. In Section III, fast and slow approximate waves in two-stream quantum 
plasmas are identified. The coupling of these waves is then analyzed in detail, in order to 
understand why some wave-numbers are stable and others are not, when varying the strength 
of quantum effects. A detailed account on the limitations of the quantum two-stream model 
for large quantum parameters H due to the neglect of collisions is included in Section IV. 
Section V is reserved to the conclusions. 



II. TIME-AVERAGED ENERGY DENSITY OF ELECTROSTATIC OSCILLA- 
TIONS 

Consider two counter propagating electron beams with equal equilibrium number densities 
n\ — Ui — no/2, and equilibrium velocities ui = —u 2 = u > 0, in presence of an immobile 
ionic background of particle density uq. Following the same notation of Ref. jlj], it is 
convenient to normalize the frequencies to the plasma frequency u)„ and the wave-numbers 
to ujp/uQ. Using the two-stream quantum hydrodynamic model a dielectric function 
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1 — F(Q) is derived, where the characteristic function F is 
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(Q + K) 2 — H 2 K 4 /A (n- K) 2 -H 2 K*/4 

in terms of a non-dimensional quantum parameter H = hw p /muQ, where h is Planck's 
constant over 2n and m is the electron mass. 

The dispersion relation e = is a second-order polynomial equation for Q 2 , with solutions 
n 2 = n%(K), where 
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corresponding to four possible branches of the eigen-frequency Q as a function of the wave- 
number K. We use parity properties to restrict the analysis to positive K and Q values. As 
detailed in Ref. when < H < 1 there is instability provided K < Ka (semiclassical 
branch) or Kb < K < Kq (quantum branch), where Ka < Kb < Kq are given by 
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On the other hand, when H > 1 the instability condition is just K < Kc (see Fig. 1 of Ref. 

In order to evaluate the energy content of each wave mode, it is necessary to consider the 
time-averaged energy density < W e > for longitudinal plasma oscillations, which is 



<W e >=j 



e d(tt£ h ) 



on 



(5) 



where Eq is the vacuum permittivity, is the Hermitian part of the dielectric function 
and E\ is the amplitude of the electric field perturbation. Since the underlying model is 
dissipation-free, one has e = Eh- Eq. (|5J) is still valid in our quantum plasma system, 
since it is derived from the Maxwell's equations only. Quantum effects are contained in the 
modified dielectric function. Moreover, < W e > includes both the contributions due to the 
electrostatic energy and the "acoustic energy", defined as the kinetic energy associated to 
the coherent particle wave motion [ll|. The expression for < W e > can equally be deduced 
from a generalized Poynting theorem, like for the classical two-stream instability [l^ . 
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Proceeding from Eqs. (pQ) and (jSJ), one get d(fl€h)/dQ ~ omitting a complicated 

positive factor, where 

^(0) = -6 Ji 4 + H 2 K 6 + + K 2 (A — H 2 K 2 ) tl 2 + 2 fi 4 . (6) 

8 

Considering f2 = Q + (K) from Eq. ([2]), one finds, 

tff{n+) = 1 + 8K 2 + 4H 2 K 6 + {l+4K 2 )Vl + 8K 2 + 4H 2 K 6 >0, (7) 

so that this mode is always a positive energy wave. 

On the other hand, taking Q = Q_ (K) from Eq. ([3]) gives, 

= 1 + 8 K 2 + 4 # 2 A' 6 - (1 + 4 K 2 ) Vl + 8 K 2 + 4 if 2 if 6 . (8) 

From Eqs. (j3J) and (jHJ), it can be shown that ^>(fl-) < if and only if K < Kc- Therefore, 
if K > Kc, the mode Q = Q_(K) is a stable positive energy mode; if K < Kc, the unstable 
modes with Kb < K < Kq (see Eq. H]) and Q = fi_ (K) are negative energy waves. Actually, 
they correspond to an absolute instability as they are of the form Q = ij, for real 7 > 0. 
The stable mode Q = tt-(K) for Ka < K < Kb, which exists only for H < 1, has negative 
energy. Finally, one has i/)(p._) = for the marginally stable wave-number K = Kc- 

Further insight can be gained analyzing the characteristic function F(Q) given by Eq. 
(CE|). It has vertical asymptotes at Q = ±f2> and Q = ±f2<, where f2> = K + H K 2 /2 and 
Q < = K — H K 2 /2. Since the dispersion relation is quadratic with real coefficients for Q 2 , 
stability is assured when the graph of F(Q) intercepts the horizontal line F = 1 four times. 
Actually, the case K = Kc is special because the quartic equation for Q degenerates into 
a quadratic one, which can be shown to correspond always to stable oscillations. Figure 1 
shows a typical unstable case when K < Ka, for K = 0.9, H = 0.1. 

On the other hand, Figure 2 explains why the wave-numbers satisfying K > Kc are 
stable, since the graph of the characteristic function always intercepts the horizontal line 
F = 1 four times. 



III. FAST AND SLOW APPROXIMATE MODES IN ELECTROSTATIC TWO- 
STREAM QUANTUM PLASMAS 

Figure 3 shows the dispersion curves for H < 1. Since the dispersion relation is an alge- 
braic equation with real coefficients, it is possible to follow Sturrock rules 0, [lo|, identifying 
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FIG. 1: Characteristic function for K < Ka- In the example, K = 0.9, H = 0.1. It correspond to 
(absolute) instability, since F(0) > 1. 



Ftfl) 




FIG. 2: Characteristic function for the stable wave-numbers K > Kq- In the example, K = 1, 
= 0.5. 



several stable/unstable zones in the if) space. Here we consider real wave- numbers and 
don't analyze the amplification problem. Curve 1 is a positive energy mode parametrized by 
fl = Q + (K) given by Eq. (j2J). Curves 2 and 3 are both described by Q = f2_(iC) given by 
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FIG. 3: Dispersion curves for H = 0.5. Curve 1 is a positive energy mode parametrized by 
Q = Q + (K). Curve 2 is a negative energy mode, while curve 3 is a positive energy mode. Both 
curves 2 and 3 are described by Q = Q-(K), see Eqs. fl2H3D. 



Eq. ([3]). However, according to the preceding analysis, curve 2 is a negative energy mode, 
while curve 3 is a positive energy mode. The coupling of these waves gives rise to the purely 
quantum (absolute) instability for large wave- numbers, Kb < K < Kc in Fig. 3. 

As apparent from Fig. (3), it is not possible to have an exact coupling between the 
negative and positive energy waves described by curves 2 and 3 respectively, since they have 
no intersection. Nevertheless, from the structure of the graphics, one might suspect that 
some approximate coupling can take place. To give support to this conjecture, a useful 
approach is to identify approximate negative and positive energy modes corresponding to 
the exact waves in an appropriate limit (e.g. for sufficiently large wave- numbers). Hence we 
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where two fast 
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approximate space-charge modes can be identified. The subscripts r and I refer to quantum 
plasma longitudinal oscillations on the rightward and leftward beams, respectively. Indeed, 
in the reference frame of the beam moving to the right, the Doppler-shifted linear terms in A 
would disappear in Eqs. ffTOj) and ffl2l) . Restoring momentarily physical units, we would then 
have u 2 = u 2 + h 2 k 4 /(4:m 2 ), which is the quantum modified Langmuir dispersion relation. 
Here, to is the wave-frequency, k is the wave-number and cj& = oj p j\[2 is the beam's plasma 
frequency. Moreover, using Eqs. (J2H3]) it can be proved that Q + ~ Qf^ r ~ K + H K 2 /2 and 
— f2_ ~ Q s>r ~ K — H K 2 /2 for large A, so that the fast and slow waves are the asymptotic 
forms of exact branches of the dispersion relation. It is in this phase velocity context that 
the (/, s) subscripts refers to "fast space-charge wave" or "slow space-charge wave", when 
focusing on a Parti cu, ar bea, Q. Si^ar considerations app, to the Leftward beam: 
0_ ~ Sl ftl ~ —K + HK 2 /2 and ~ tt Stl ~ -A' -HK 2 /2 for large A. 

Calculating the time-averaged energy density < H4 > for electrostatic oscillations using 
Eq. (J3J), it is directly verified that VLf r and Q s j are positive energy modes for all K. On the 
other hand, Qfj and f2 Sjr can have negative energy content depending on special conditions. 
As for the classical two-stream instability [3], instability is expected when a wave on a 
beam has free energy to drive an instability in the counter propagating beam. Inspecting 
Eqs. ( TTTJHT3T) shows that the only possible couplings are ttfj = Q Str or, reversing directions, 
Qf t1 . = Q s j. Focusing on A > and positive frequencies, instability is expected when the 
fast positive energy wave of the rightward beam couples to a negative energy wave on the 
leftward beam. It is found that the slow wave Q Sjr is a negative energy mode provided 
A < Ac (assuming non-negative £l 3)r )- Hence the matching condition is Q,fj = Q s ,r, or 
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as illustrated in Fig. 4, which shows intersection of the fast leftward and slow rightward 
modes for H = 0.1. As in the classical case, the coupling occurs for Re(fi) = 0, in analogy 
with ideal MHD instabilities [l2]. 

In the context of this interpretation, the wave-numbers K m satisfying the coupling con- 
dition ( Ti~4l) correspond to maximal instability growth-rate. For H ^ 0, they are given by 
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Notice that in the classical case where H = 0, Eq. ([14]) yields only the solution K 
Taking the plus sign in Eq. ffTS]) . one finds 
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Assuming H 1, we get to leading order, 



j + 0(H 3 ) >K B ^K c -^ + 0(H 3 ) 



(17) 



where Kb and Kc are defined in Eq. (jlj). Therefore, Kb < K m ^ q < Kc, which is exactly 
what should be expected for an instability arising from the coupling of the positive energy 
mode shown in curve 3 and the negative energy mode shown in curve 2 of Figure 3. This 



explains the physical origin of the purely quantum instability described in Ref. occurring 
for large wave- numbers, when H < 1. For a fixed value of H < 1, there is a good agreement 
between K m ^ and the exact wave-number for fastest growing quantum instability. The 
discrepancy comes from the fact that Eqs. ffTUHTB"]) show just approximate modes. 
On the other hand, taking the minus sign in Eq. (TTo*j) gives 

1/2 

(18) 
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We get, to leading order, 



K m , c = + ^= + O(H') <K A ^l + ^ + 0(H*) , (19) 

where is defined in Eq. fl4]). The wave-number K m c can be interpreted as the semiclas- 
sical branch, since it corresponds to the exact classical wave-number for maximal instability, 
K c = l/v2- Moreover K myC < Ka corresponds to the coupling of the positive (curve 1) and 
negative energy (curve 2) branches in Figure 3. Finally, we found a satisfactory agreement 
between the exact and approximate values of the wave-number for maximal growth-rate, at 
a fixed H. 

When H > 1, the elliptic-like branch of Figure 3 disappears. The instability for K < Kc 
and H > 1 can also be understood in terms of negative energy modes. However, we omit 
these considerations since these parameter ranges are outside the scope of the quantum two- 
stream model. Indeed, as shown in the next Section, H should be small to avoid the strong 
coupling regime. 



IV. NECESSARY VALIDITY CONDITION FOR THE COLD QUANTUM TWO- 
STREAM MODEL 

In a cold two-stream plasma model, be classical or quantum, it is implicit that each beam 
has a well-defined identity, which is reasonable only if the respective velocity dispersions are 
much smaller than the beam velocities. It is relevant to analyze the implications of this 
point in view of the collisionless hypothesis underlying the original equations [l|. 

Consider first the case of classical statistics, which applies for sufficiently small particle 
densities, or large temperatures. In this case, the beam velocity dispersion is the thermal 
velocity v t h, so that the present model assumes u ^> u th , which is equivalent to a cold 
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beam approximation. Moreover, the thermal velocity of the beams is a free parameter, 

1 /3 

independent of the particle density. It can be verified that gc = e 2 n /(eomv 2 h ) <C 1 
(small classical coupling constant), together with H > 1, implies some relativistic velocities 
which are out of the scope of the model. Therefore, H -C 1 is necessary to support the 
collisionless hypothesis. 

Let us now rule out the possibility of a large quantum parameters for dense degenerate 
plasmas a well. In this case, due to Pauli blocking, the system is all the more ideal that the 
density is large, and it is appropriate 115] to define the quantum coupling parameter gq = 
e 2 n^ 3 1 (eq Tp) ~ {hujp/ (mu 2 F )) 2 . Here, Tp and uf are the electronic Fermi temperature and 
velocity, respectively. In the degenerate case, up measures the velocity spread of each beam, 
and increases with the density. Hence, the quantum two-stream model assumes u ^> up, 
imposing a bound on the electron particle densities. However, for such small densities 
one has H <C hu p /(mu 2 F ) ~ ^ <C 1. Here again, the collisionless and the well-defined 
two-stream hypothesis, demand H -C 1. The present work provide physical arguments to 
understand the quantum instability branch for small H and large wave-numbers, where the 
weak coupling regime apply. In addition, large quantum parameters are admissible in the 
context of a quantum Dawson iV— stream model , in which case the beam velocities 

need not to be large in comparison to the Fermi velocities. 



V. CONCLUSION 



The main results of this work stem from the factorized form of the dispersion relation, 
as expressed in Eq. ([9]). This allows to identify fast and slow waves propagating in both 
the positive and the negative directions. The quantum two-stream instability grows when 
the free energy available in a positive energy wave carried by one beam, is transferred to 
a negative energy wave carried by the other beam. The coupling of these waves and their 
energy exchange has been discussed in terms of their electrostatic energy density, giving 
rise to stable or unstable scenarios. While the mathematical techniques for the quantum 
plasma problem are the same as for the classical two-stream plasma, the quantum dispersion 
relation is more subtle as evidenced by Fig. ([3]). The above analysis can in principle be 
pursued on similar problems such as the quantum beam-plasma instability and on theories 
with a discrete structure such as the quantum multi-stream model. 



11 



Acknowledgments 

This research was financially supported by the Alexander von Humboldt Foundation, by 
projects FIS 2006-05389 of the Spanish Ministerio de Educacion y Ciencia and PAI-05- 
045 of the Consejeria de Educacion y Ciencia de la Junta de Comunidades de Castilla-La 
Mancha. We acknowledge the anonymous referee for his constructive questions, in particular 
for pointing out the model's applicability condition on the parameter H. 



[1] F. Haas, G. Manfredi and M. Feix, Phys. Rev. E 62, 2763 (2000). 

[2] D. Anderson, B. Hall, M. Lisak and M. Marklund, Phys. Rev. E 65, 046417 (2002). 

[3] F. Haas, G. Manfredi and J. Goedert, Braz. J. Phys. 33, 128 (2003). 

[4] S. Ali and P. K. Shukla, Eur. Phys. J. D 41, 319 (2007). 

[5] A. Mushtaq and R. Khan, Phys. Scr. 78, 015501 (2008). 

[6] H. J. Ren, Z. W. Wu, J. T. Cao and P. K. Chu, Phys. Lett. A 372, 2676 (2008). 

[7] A. Bret, Phys. Plasmas 15, 022109 (2008). 

[8] A. Bret, Phys. Plasmas 14, 084503 (2007). 

[9] P. A. Sturrock, Phys. Rev. 117, 1426 (1960). 
[10] A. I. Akhiezer, I. A. Akhiezer, R. V. Polovin, A. G. Sitenko and K. N. Stepanov, Plasma 

Electrodynamics (vol. 1, Linear Theory), Pergamon Press, Oxford, 1975. 
[11] T. H. Stix, Waves in Plasmas, Springer Verlag, New York, 1992. 
[12] C. N. Lashmore-Davies, Phys. of Plasmas 14, 092101 (2007). 

[13] R. J. Briggs, in Advances in Plasma Physics (vol. 4), ed. A. Simon and W. B. Thompson, 

John Wiley & Sons, New York, 1971. 
[14] L. D. Landau and E. M. Lifshitz, Statistical Physics, part 1, Butterworth-Heinemann, Oxford, 

1980. 

[15] See Eq. (2.12) of G. Manfredi, Fields Inst. Commun. 46, 263 (2005). 
[16] J. Dawson, Phys. Fluids 4, 869 (1961). 



12 



